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Abstract 

The ongoing and planned experimental activities with direct reference to light unflavoured pseudoscalar mesons 
motivate a new theoretical study regarding their properties. An overview including details on new precise calculations 
is presented. 
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1. Introduction 

The subjects of this work, light pseudoscalar mesons, play a prominent role in hadronic physics. Our even more 
focused interest into the unflavoured particles, namely n°, r\ and eventually 77', is motivated by the wish to avoid 
a discussion on K° decays. Of course, not because they are not interesting, but such decays violate hypercharge 
conservation and are suppressed by Fermi coupling constant G 2 F , meaning that studying only unflavoured ones enables 
to reduce standard model to QCD and simplifies the problem. For studying QCD at low energy region, in our case 
enlarged at most only by QED corrections, a standard tool successfully developed in recent years is called chiral 
perturbation theory (ChPT). In this contribution we will mainly discuss properties of n° as many radiative 77 decays 
are technically very similar and what one obtains for n° can be simply converted also for 77 decay prediction. On the 
other hand, jt° being the lightest meson cannot decay into other hadronic states, therefore a hadronic discussion for 77 
decays is inevitable. As an important example of such processes we will briefly mention 77 — > 3n decays. 

The decay modes of 7r° were subjects of many experiments in the past (including e.g. SINDRUM coll. at PSI), 
present (e.g. KTeV, or PrimEx at JLab) or future (NA62 at CERN). Using the conserved vector current hypothesis 
we can connect the vector form factor (i.e. charged pion) to the lifetime of the neutral pion (cf. PIBETA [ 1 1). Exper- 
iments have reached (or plan to reach) a level of precision which makes it mandatory to reopen previous theoretical 
calculations to achieve appropriate order (NLO or NNLO). This can, on one hand, help us to verify and fix the un- 
derlying structure of the low energy effective theory of QCD - ChPT (e.g. pion decay constant, low energy constants, 
power-counting, etc.). On the other hand it can set a framework for studying new physics beyond the SM (e.g. KTeV's 
discrepancy with a theory for n° — > e + e~). 
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As mentioned the 77 meson can be treated technically very similarly. However, due to its mass one can also study 
different leptonic variants and combination (for example // in place of <?) and last but not least its hadronic decays. 
This can provide us with important information on isospin breaking effects and again test the internal consistency of 
ChPT. 

We will mainly focus on four most important allowed decay modes of 7r°: yy, e + e~y, e + e~e + e~, e + e~ (with branch- 
ing ratios 0: 0.98798(32), 0.01 198(32), 3.14(30) x 10~ 5 , 6.46(33) x 10~ 8 , respectively). For this purpose one can use 
two-flavour chiral perturbation theory (ChPT, for a review see e.g. [3 |) which can simply incorporate corrections to the 
current algebra result attributed either to m u j masses or electromagnetic corrections with other effects hidden in the 
low energy constants (LECs), denoted by at next-to-leading order (NLO). However, phenomenologically richer 
S U(3) ChPT must be also employed in order to obtain a numerical prediction. This is especially true for the studied 
anomalous processes as in this case the initial symmetry for the two flavour case must be extended and the number of 
monomials in SU(2) increases |4|. 

We will not limit our focus only on "standard on-shell" decays. It is clear that both on-shell and off-shell or semi- 
of-shell vertices, especially 7r ow y*y ( *\ play a crucial role in many other experiments, from the famous g — 2 (cf. 0) 
via virtual photons stemming from e + e~ (see e.g. the recent paper |6|) to astrophysics. Our first aim is the common 
formulation of these interrelated processes in the given formalism at the given order (either NLO or NNLO) motivated 
by the precision of the appropriate experiments. 

2. n -* yy 

The 7T° meson has a prominent position among all hadron particles as being the lightest state of them. Its primary 
decay mode is thus n° — > yy which is connected with the famous Adler-Bell-Jackiw triangle anomaly [7|. It saturates 
the decay width with almost 99% and plays an important role in the further decay modes (see the following sections). 
The history of 7r° — > yy is going back to Steinberger's calculation JSJ, for a review see e.g. (9|, and also [ 101 . The 
prediction estimated from the chiral anomaly using current algebra agrees surprisingly very well with experiments. A 
first attempt to explain the small existing deviation from the measurement was made by Y. Kitazawa ifTTI . At that time 
a new experimental prediction from CERN-NA030 iTPUl suggested a smaller value for the partial width 7.25 ± 0.23 
eV (statistical and systematic errors combined in quadrature). Older experiments (Tomsk, Desy and Cornell |13|), 
seemed not to be so precise (7.23 ±0.55, 1 1.7+ 1.2, 7.92 ±0.42 eV, respectively); they relied on the so-called Primakoff 
effect [ 14 1 that is based on measuring the cross section for the photoproduction of the meson in the Coulomb field. The 
more precise number from the direct measurement at CERN motivated Y. Kitazawa to explain the 8 ~ 9% discrepancy 
by including QED correction and the 77/77' contribution. These corrections were not, however, large enough to explain 
the discrepancy which was attributed by the author to a possible 7r(1300) contribution. Furthermore, it was found out 
in this work that the contribution from multi-pion states must be small. This was verified explicitly also within ChPT 
with the remarkable observation |[T5l that at one-loop order there are no chiral logarithms (either from pions or kaons). 
The n-rj-r]' mixing and electromagnetic correction were reconsidered relatively recently in ifTBI . 

The spread in the data basis of the PDG, summarized in the previous paragraph, shows, however, that the quoted 
errors seem to be underestimated [17|. The present situation fortunately looks more optimistic as the world average 
accuracy of 8 % is planned to be improved to the level of one or two percents in ongoing experiment PrimEx at JLab. 
The official Run-I result quoted in [18] is r(7T° — > yy) = 7.82 eV ± 2.8%. This was the main motivation for a new 
study of 7T° — > yy in 1 19 1. The correction to the chiral anomaly due to the finite mass of light quarks was reconsidered 
using strict two-flavour ChPT at NNLO. We will summarize here this remarkably simple result (note that it involves 
a two-loop calculation and that it represents formally a full 0(p 8 ) result). Defining a reduced T amplitude 

A = e l s llva f S efe™]i[k v 1 T , (1) 



we have for the partial decay width 



r yr = ^W \t\ 2 



(2) 
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Up to and including next-to-next-to-leading order corrections 



F„T NNW = ^ + ^ (-<'' - 4cf + cY{) + ^HU - rn n H5< f : «•» ' : 2< f ) 
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(3) 



where the chiral logarithm is denoted by L„ = log and A + , A-, A can be expressed as follows in terms of 

renormalized chiral coupling constants (d Wr refer to combinations of couplings from the NNLO Lagrangian, i.e. of 
order p 8 in the anomalous sector), 



1 



-dT Ou) - 8c' 



HQ 2 



5127T 4 
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"288 



16 
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A_ = ^ [^00 + ^ (5cf + cf + 2cf0] 



= dfty) - 128F 2 / 7 (c 3 H "' + c7 r ) , 



(4) 



with Riemann zeta and Clausen function: £(3) = 1.202.. and Cl2(?r/3) = 1.014.., respectively. 

All effects that were carefully studied e.g. in ifTTIl and ITS! are now hidden in the LECs and chiral logarithms L n 
(with the exception of QED corrections that must be added by hand to the latter formula, see Ifl9l ). For a detailed 
phenomenological study of the previous formula see lfT9ll , our best estimate has led to 



yy 



(8.09 ± 0.1 l)eV. 



(5) 



The importance of various input parameters is depicted in Fig. [T] One can see immediately the importance of F n in 




change in [%] 



Figure 1: The dependence of n° — > yy decay width on various parameters, namely S 1/(3) LEC CF, R and F n (see main text). 

n — > yy decay. The F n , on the other hand, is determined from the weak decay of n + based on the standard V- A inter- 
action. The new proposed variant of this interaction beyond SM assumes contributions of right-handed current which 
would lead to a change of F„ [20|. Determination of this constant directly from n° lifetime can provide constraints 
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on such contributions. This again put big effort to minimize the uncertainty stemming from other parameters as for 
example visualized in Fig.jl] One of them, a quark mass ratio R will be subject of Sec. [7] the second one reflects 
the intrinsic connection with 77 — > yy decay as calculated in three-flavour ChPT. However, better understanding of 
77 — > yy again in full two-loop calculation is important and probably necessary as motivated and explained in ED . 



3. n -* e + e y 

If 7t° — > yy represented 99% of all decay modes, 7r° — > e + e~y represents again more than 99% of the rest modes. 
It is thus the second most important decay mode with a branching ratio ~ 1.174 ± 0.035% and nowadays called Dalitz 
decay, after R.H. Dalitz who first realized its connection with two-photon production [22 1 . Knowing the branching 
ratio one can also use, at least in principle, the Dalitz decay to extrapolate the total decay width, which can serve as an 
independent possibility how to measure the life-time of tt° (however, handicapped by the larger dependence on error 
of this ratio). This year's change of the official PDG number is due to ALEPH archived data on hadronic Z boson 
decay which has enabled to reconstruct 12,490 Dalitz decays with a result 

ALEPH: r^^/Tyy = (1.140 ± 0.041)%, (6) 

that led to the update 

PDG 10: r^^/Tyy = (1.188 ± 0.035)%. (7) 
To conclude let us also summarize a theoretical prediction: 

theory: T^^^/T^ = (1.1851 + 0.0104 + 0.0018)% = (1.1973 + 0.0055)% , (8) 

where the first number stands for the leading order, predicted already by Dalitz [22], the second represents the radiative 
corrections (numerically first done in [23 1) and the last number [24] stands for the prediction of QCD corrections 
and two-photon exchange contribution (which was neglected previously in QED corrections). The final error was 
estimated as a half of all QED corrections. 

The Dalitz decay is very often used experimentally as the normalization mode not only for rare pion modes 
(see also below) but also for kaon decay modes, and thus its precise value has impact on these measurements. Its 
uncertainty has in fact a direct effect on external systematic error and different central values can produce substantial 
shifts in the final predictions. 

The motivation for studying the decay width of 7T° — > e + e y is thus two-fold. The precise and well-understood 
theoretical prediction with model-dependent QCD contribution suppressed by phase-space integration can serve as a 
calibration of the experiment. On the other hand, if theory is well under control and still some measurement would 
signalize some discrepancy the theory is missing something new and important. 

The total decay rate is, however, not the whole story. It turns out that the corrections to the differential decay 
which were taken as negligible are indeed important. The reason is that there is a part of the phase space where, 
roughly-speaking, the correction to the differential decay width is positive and a part where it is negative; and only 
summing these parts together gives us the small number in (|8j. It is clear now, that in physically relevant applications, 
when we have to cut some parts of the phase space, these corrections can become important. Let us discuss this a little 
bit more in detail (for details see 11241 ~). First it is useful to define two kinematic variables that represent the normalized 
di-lepton invariant mass and the difference of energy for positron and electron (normalized to the photon energy in the 
pion rest frame): 

(P+ + P-) 2 2 „ „ , 2 4 ™ 2 

X = ; , V < X < 1 , V = r-, 



IP ■ (p+ - »_) , , / 4m 2 

2 „ n o ^ 
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(m denotes the electron mass, p 2 — p 2 — m 2 and P stands for the pion momentum). The next-to-leading corrections 
to the differential decay rates can be described as 

dr dr L0 

5(x,y) ■ 



dx dx 

where the corresponding LO partial decay rates have a relatively simple form 

dr LO _ a 3 M ff o (1 -x) 3 2 2 



dxdy dxdy ._. 

dr dr L0 



dxdy (4tt) 4 F% x 1 " 

cr e (xM 2 )(xM 2 +2m z ) (10) 



dT L0 _ a 3 SM^o (1-x) 3 f - w , i o 2 



dx (4tt) 4 3 F 2 x 2 



(in fact, integrating the last equation one can verify the Dalitz result, i.e. the first number in ([H}). With these quantities 
in hand we can extract information on the QCD part of the form factor TjPyyiq 2 ), which is related to the doubly 
off-shell Tiyy transition form factor defined as 



^ d«xe Ux (0\T(f(x)f(0)\7i v (P)) = -ieT^PpT^y^AP - If), (11) 
by 

T^yyiq 2 ) = T„o rr (0,q 2 ). 

The Dalitz decay can provides us with the information on the transition form factor in the time-like region. This is 
usually specified by its slope parameter a„ 



Tfyfisft = ^o yy ,(0) [I + a n ^ +■■■}, (12) 

1V1 _ n 



and from the experiment can be obtained by subtracting the QED corrections via 

dr"'' dr L0 dr L0 



- Sqed(x) — — = — — [1 +2xa„], (13) 



dx dx dx 

The direct measurements in the time-like region is endowed with large errors 

SaclayED: fljr = -0.11 + 0.03 ± 0.08 
TRIUMF J26): a„ = +0.026 ± 0.024 ± 0.0048 
PSI|27): a n = +0.025 + 0.014 ± 0.026, 

whereas the values extracted from the extrapolation of data at higher energies in the space-like region, Q 2 = -q 2 > 0.5 
GeV 2 are more precise: 



CELLO El: a„ = +0.0326 + 0.0026 + 0.0026 
CLEO 11291 : a„ = +0.0303 + 0.0008 + 0.0009 + 0.0012 



The theoretical calculation l24l (here without isospin correction for simplicity) is given by 



32tt 2 M 2 w , M\ . M 2 1 a 

n -c^(M 2 v ) f- l+21n^ - — - (14) 

3 13 V vJ 96n 2 F^ M 2 V > 360 tt V 
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We can see that the EM corrections are very smalj^J which signalized its tight connection with a low energy constant 
(or via S 1/(2) - S U(3) relations IE) to S i/(3) LEC C™. The LMD prediction for eg leads to 

theory (24]]: a n = +0.029 ± 0.005 (15) 

Let us also note that the off-shell form factor f^y plays an important role in the 7r°-exchange contribution to a 
hadronic light-by-light scattering contribution in g - 2 and thus more experimental information on this factor can help 
us to understand the consistency in different approaches. 

4. 7i° — » e+e~e + e~ 

Having the well established decay of tt° into two photons it is clear that the pion cannot be a J — 1 state (so- 
called Young-Landau theorem ll30l ). However, using directly n° — » yy to verify whether it is a (pseudo)scalar is 
experimentally impossible. It was thus suggested in l3T| to use the double-internal conversion, the so-called double- 
Dalitz decay. The experiment was performed at Nevis Lab 1321 in a bubble chamber with the following result for the 
branching ratio: 

pPDG'08 

-^^1 = (3.14 + 0.30) x 10~ 5 (16) 

1 tot 

which was used as only relevant measurement for almost half century. The experiment also confirmed the negative 
parity of n° known from the previous indirect studies via the cross-section of tT capture on deuterons. However, the 
significance of this direct measurement was only 3.6 cr. Recently the long standing experimental gap was filled with 
a new measurement in the KTeV-E799 experiment at Fermilab [33] giving a branching ratio (including the radiative 
final states above a certain cut as tacitly assumed for all Dalitz modes) 

pKTeV 

-^p!£l = (3.4 6 + 0.19) x 10~ 5 , (17) 

1 tot 

which is in good agreement with the previous experiment. In addition to the precisely verified parity of n° (which 
represents its best direct determination) this experiment sets the first limits on the parity and CPT violation for this 
decay. More precisely, having a n°y*y* vertex C IJ y po -F f ' y F pa 'n we can study, using the following decomposition (for 
details see Il3"4"l ) 

+ sin ^(g^gya- - g^o-gyp) , 

the parameters £ and 6 which represent parity mixing and CPT violation parameters. For details see 0311 : for example 
their limit on the mixing assuming CPT conservation is ( < 1.9°. 

A detailed analysis of the radiative corrections in |34] showed that they seem to be very important in extracting 
physically relevant quantities. This motivates us to reopen this subject l35l in the same manner as was done in [24|. 
The simply looking task of attaching another Dalitz pair on the virtual photon line is complicated (in the defined power- 
counting) by the necessity to include a pentagonal diagram [34]. This strengthens the need of a correct description of 
the off-shell ji°y*y* vertex, which can be, on the other hand, directly studied in the next mode. 

5. n° e + e~ 

In the previous decay modes the fully off-shell n°y*y* vertex was suppressed by the dominant semi-on-shell 
contributions. As this is not true anymore for 7r° — > e + e~ it naturally represents the simplest and cleanest candidate for 
studying not-well understood effects of QCD, i.e. T^yyik 2 + 0, 1 2 + 0), even though the process itself is suppressed 
by approximate helicity conservation and two extra powers of a. This is supported by the existing experiment at 



1 Let us note that in previous experimental analyses the two-photon contribution were neglected and their omission would lead to approximately 
0.005 correction into the right direction towards the independent CELLO or CLEO result. 
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Fermilab (KTeV E799-II) 11361 . Comparing with the previous measurements their result has increased significantly 
the precision and provide thus the most important contribution to the present PDG's average 

pPDG 

-p^- = (6.46 ± 0.33) x 1(T 8 . 

1 tot 

This process was first calculated in 11371 and proceeds, as mentioned, via two intermediate photons and at LO is thus 
represented by an one-loop (triangle) diagram. One can get (for details see e.g. (38]) 

^=2(^)V(M^(M2 )| 2 , (18) 

lyy X 7TM n / 

with the amplitude given by 

3l( S )=X(M)-l + l^+C(s), (19) 
2 2 

where C represents the scalar one-loop triangle. The imaginary part can be calculated in a model independent way 
by cutting the photon lines and knowing the on-shell form factor T n oy, y ,(Q, 0) (in fact as we have a normalization to 
two-photon decays it has to be equal 1) with the result 

7T 1 — cr 

Imyi= — In , (20) 

2cr 1 + cr 

implying thus a unitarity bound on the ratio r e + e -/Tyy > 5 x 10~ 8 . The real part can be now calculated via a dispersive 
integral leaving us with one unknown subtraction, or equivalently by x(m)- The techniques of large 7Y C together with 
the LMD approximation (V represents the p meson) leads to [38 1 

j | j?2 pth.no-rad 

X(fi = M v )= —-4n 2 -^ =2.2 + 0.9 => -^f — = (6.2 ± 0.3) x 10~ 8 . (21) 
This should be compared with the experiment after removing the effects of final state radiation: 

pKTeV.no-rad 

= (7.57 + 0.39) x 10~ 8 , (22) 

I yy 

which is 3.5cr off from the mentioned theoretical prediction. 

At present the theoretical activities concerning this process focus on two main directions: z) understanding the 
discrepancy within the SM - i.e. calculating radiative corrections or employing a different model for the QCD part 
1391 . ii) It naturally represents an ideal candidate for testing new models beyond SM, it can set valuable limit on light 
dark matter scenarios, supersymmetric extensions (axion), etc... 

Before concluding let us also mention the weak contribution. This process has in fact a direct tree level contribution 
mediated via n° — > Z* — » e + e~ . The amplitude would be proportional to the Fermi constant: 



^Jteu- ~G F F K muy 5 v, (23) 



which makes it three orders of magnitude smaller than the dominant EM contribution. Beyond-SM scenarios would 
use a similar relation (with different coupling Gf ) for introducing the effect of a light vector particle (e.g. U-boson). 



6. n -* vv, invisible, extra-light particles 

As studied for example by |40| there is a tight connection for 7r° — » vv with cosmology so that the strong limit on 
this decay obtained are much higher than those in the laboratory: 2.7 x 10~ 7 (E949 based on K + — > n + n° fiTTl ). This 
decay mode represents not only two neutrino decay modes but all possible combination and weakly coupled exotics 
and also even more generally tt — > invisible. 
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Within the SM (extended by massive neutrinos) one can use the calculation of the weak sector for electron from 
previous section with a prediction 

•HjSlUr- = ^2G F F„myuy 5 v, (24) 
with one subtlety - we don't know the mass of the neutrino and its nature (for the Majorana type the amplitude is 
twice bigger). The relative branching ratio is then (in the Dirac case) 

T V y ( 4nFlG F m v s2 f 4mf 

rj^ = ' a ~M' V Af* (25) 

With a direct limit on the absolute tau neutrino mass m T v < 18.2 MeV we get a reasonably high limit on branching 
ratio < 5 x 10~ 10 (and twice bigger for Majorana case). Note the maximum for the ratio for m v = M/ Vf>, which is, 
however, ruled out. 

The helicity suppression for ji° — > vv can be avoided in the decay mode 7r° — > vvy. This decay mode is also 
interesting as it depends on the actual number of light neutrinos. 

A general possibility to set the constraints on an extra-light long-lived (/non-interacting/weakly interacting) neutral 
vector particle X via decays of known particles opens naturally this question also for the exotic 7T° decay. 



7. 17 -> 3tt 

We have selected this representative of 77 decay modes as the most important example of the hadronic decays 
for the studied particles. A brief look into the history of this mode ([42 1 for NLO and [43 1 for NNLO) shows 
even after two-loop calculation a discrepancy between theory and experiments. This could be attributed either to 
bad convergence of the ChPT series and necessity to somehow (at least partially) re-sum higher orders or to some 
problems in the calculation at NNLO (e.g. to a wrong estimate of LECs). The new ongoing efforts represent in some 
sense combinations of these two possibilities ( [ 44 1 , [45|, |46|). 

Apart from the obvious motivation to better understand how to perform NNLO calculation in ChPT in the three 
flavour case, which seems to be a problematic subject we have in hands a process which vanishes in the isospin 
limit. It means that a precise measurement together with a good understanding of the theory should provide us with 
parameters which can quantify the isospin breaking, for example: 

m s -in r, , / ra, >. 

R=— Q 2 = Ul + ^]R. 

rrid — m u y m ' 

The absolute value of the partial decay width for 77 — » 37r is experimentally obtained via normalization to 77 — > yy. 
A change in one decay width has thus influence in other (a change by 1% in T(rj — > yy) input shifts R by w 0.2). 



8. Summary 

The new experimental activities in the low energy physics that concern directly n° or 77 decay modes call for a 
more detailed theoretical study in this area. We have discussed some allowed decay modes that represent important 
tools for studying basic phenomena of the underlying theory: QCD. Namely, jt° — > yy and 77 — > yy played an 
important role in understanding a symmetry pattern of the theory as they are directly connected with the so-called 
U(1)a anomaly. QCD enlarged by photons possesses, however, two such anomalies. The first one, internal, connected 
with QCD only, proportional to gluonic term G /JV G /JV , dubbed [/(l)-problem and the resulting strong CP problem is 
still an open issue. As a remnant of the anomaly, the 77' plays a more important role than naively expected and has to 
be included in a theoretical consideration. The second anomaly, external, in our case connected with electromagnetic 
interaction (or F^F^) explains why n — > yy can decay so quickly even it should be suppressed due to Sutherland's 
theorem. Furthermore, 77 — > yy and 77 — > 3n represent 95% of all 77-decay modes and are thus perfectly suited to study 
directly properties of the 77. Simultaneous treatment of two-photon 7r° and 77 decays, apart from testing or fixing our 
understanding of 77', can provide valuable information on the decay constants F„ and F n or quark mass ratio. 

A common treatment is thus useful and important in order to understand all phenomena. In this work we have 
mainly focused on chiral and QED corrections in order to prepare the ground for the discussion of non-perturbative 
effects or eventual new physics. 
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